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Abstract 

We consider the prompt photon production at high energy hadron colliders in the frame- 
work of fey-factorization approach. The unintegrated quark and gluon distributions in a 
proton are determined using the Kimber-Martin-Ryskin prescription. The conservative er- 
ror analisys is performed. We investigate both inclusive prompt photon and prompt photon 
and associated muon production rates. In Standard Model such events come mainly due 
to Compton scattering process where the final heavy (charm or bottom) quark produces a 
muon. The theoretical results are compared with recent experimental data taken by the 
D0 and CDF collaborations at Fermilab Tevatron. Our analysis also covers the azimuthal 
correlations between produced prompt photon and muon which can provide an important 
information about non-collinear parton evolution in a proton. Finally, we extrapolate the 
theoretical predictions to CERN LHC energies. 



1 Introduction 

It is well known that production of prompt (or direct) photons at high energies has 
provided a direct probe of the hard subprocess dynamics, since produced photons are largely 
insensitive to the final-state hadronization effects. Hadroproduction of prompt photons has 
been studied in a number of experiments [1-5]. Usually photons are called "prompt" if 
they are coupled to the interacting quarks. In the framework of Quantum Chromodynamics 
(QCD) the dominant production mechanism for the prompt photons at Tevatron and LHC 
colliders is the Compton scattering gq — > '-/q [6]. It is clear that cross section of such 
processes is sensitive to the gluon distributions in a proton. Also observed final state photon 
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may arise from so called fragmentation processes [7], where a quark or gluon is transformed 
into 7. The cross sections of such processes involve relative poorly known parton-to-photon 
fragmentation functions [8]. However, the isolation criterion which is usually introduced in 
experimental analyses substantially reduces the fragmentation component (see, for example, 
Ref. [9]). Therefore for the theoretical description of prompt photon production at Tevatron 
the detailed knowledge of parton (quark and gluon) distributions in a proton is necessary. 

Usually quark and gluon densities in a proton are described by the Dokshitzer-Gribov- 
Lipatov-Altarelli-Parizi (DGLAP) evolution equation [10] where large logarithmic terms pro- 
portional to In /i 2 are taken into account only. The cross sections can be rewritten in terms 
of process-dependent hard matrix elements convoluted with universal quark and/or gluon 
density functions. In this way the dominant contributions come from diagrams where parton 
emissions in initial state are strongly ordered in virtuality. This is called collinear factor- 
ization, as the strong ordering means that the virtuality of the parton entering the hard 
scattering matrix elements can be neglected compared to the large scale \i. 

However, at high energies (or small x ~ f^ 2 /s <C 1) effects of finite virtualities and 
transverse momenta of the incoming partons may become more and more important. These 
effects can be systematically accounted for in the /c^-factorization QCD approach [10-14]. 
Just as for DGLAP, in this way it is possible to factorize an observable into a convolution of 
process-dependent hard matrix elements with universal parton distributions. But as the vir- 
tualities (and transverse momenta) of the emitted partons are no longer ordered, the matrix 
elements have to be taken off-shell and the convolution made also over transverse momentum 
kx with the unintegrated (i.e. /^-dependent) parton distributions. The unintegrated parton 
distribution f a (x, k 2 -) determines the probability to find a type a parton carrying the longi- 
tudinal momentum fraction x and the transverse momentum ky. In particular, usage of the 
unintegrated parton densities have the advantage that it takes into account true kinematics 
of the process under consideration even at leading order. 

The unintegrated parton distributions f a ( subject of intensive studies [15, 16]. 

Various approaches to investigate these quantities has been proposed. It is believed that at 
assymptotically large energies (or very small x) the theoretically correct description is given 
by the Balitsky-Fadin-Kuraev-Lipatov (BFKL) evolution equation [17] where large terms 
proportional to Inl/x are taken into account. Another approach, valid for both small and 
large x, have been developed by Ciafaloni, Catani, Fiorani and Marchesini, and is known as 
the CCFM model [18]. It introduces angular ordering of emissions to correctly treat gluon 
coherence effects. In the limit of asymptotic energies, it almost equivalent to the BFKL [19- 
21], but also similar to the DGLAP evolution for large x ~ 1. The resulting unintegrated 
gluon distribution depends on two scales, the additional scale q is a variable related to the 
maximum angle allowed in the emission and plays the role of the evolution scale [i in the 
collinear parton densities. 

The two-scale involved unintegrated parton distributions it is possible to obtain also from 
the conventional ones using the Kimber-Martin-Ryskin (KMR) prescription [22]. In this 
way the \x dependence in the unintegrated parton distribution enters only in last step of the 
evolution, and single scale evolution equations can be used up to this step. Such procedure 
can be applied to a proton as well as photon and is expected to account for the main part 
of the conventional next-to-leading logarithmic QCD corrections. The KMR-constructed 
parton densities were used, in particular, to describe the heavy quark production in 77 
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collisions at CERN LEP2 [23], prompt photon photoproduction at DESY HERA [24] (both 
inclusive and associated with hadronic jet) and inclusive prompt photon hadroproduction at 
Fermilab Tevatron [25]. 

Recently new experimental data on the p + p — * 7 + X process at Tevatron were obtained 
by the D0 [1, 2] and CDF [3, 4] collaborations. It was found [1-4] that the shape of the 
measured cross sections as a function of photon transverse energy E T is poorly described 
by next-to-leading order (NLO) QCD calculations: the observed E T distribution is steeper 
than the predictions of perturbative QCD. These shape differences lead to a significant 
disagreement in the ratio of cross sections calculated at different center-of-mass energies 
^Js = 630 GeV and y/s = 1800 GeV as a function of scaling variable xt = 2E T /y/s. The 
disagreement in the xt ratio is difficult to explain with conventional theoretical uncertainties 
connected with scale dependence and parametrizations of the parton distributions [2, 3]. 
However, such discrepancy can be reduced [26, 27] by introducing some additional intrinsic 
transverse momentum kr of the incoming partons, which is usually assumed to have a 
Gaussian-like distribution [26]. The average value of this kx increases from (kr) ~ 0.5 
GeV to more than (kx) ~ 2 GeV as the \/s increases from UA6 to Tevatron energies [27]. 
However, in this case the full kinematics of the subprocess is not taken into account, as it 
was argued in Ref. [25]. 

The inclusive prompt photon hadroproduction at Tevatron in the /c T -factorization QCD 
approach was considered in Ref. [25] . The unintegrated parton distributions in a proton were 
obtained using the KMR formalism. The role of the both non-perturbative and perturbative 
components of partonic transverse momentum kx in describing of the observed E T spectrum 
was investigated. However, the KMR unintegrated parton densities were obtained in the 
double leading logarithmic approximation (DLLA) only. Also in these calculations the usual 
on-shell matrix elements of hard partonic subprocesses were evaluated with precise off-shell 
kinematics. 

In the present paper we apply the KMR method to obtain the unintegrated quark and 
gluon distributions in a proton f a ( ) independently from other authors. Then, we 

study inclusive prompt photon hadroproduction at Tevatron in more detail. We calculate 
the double differential cross sections dajdE^drp at two different center-of-mass energies 
\fs = 630 GeV and y/s = 1800 GeV and compare our theoretical results with the recent D0 
and CDF experimental data [1-4]. In order to estimate the theoretical uncertainties of our 
predictions we study the renormalization and factorization scale dependences of calculated 
cross sections. Also we study the ratio of cross sections calculated at different center-of-mass 
energies y/s = 630 GeV and y/s = 1800 GeV. After that we extrapolate our predictions 
to LHC energies. In all these calculations we use the expressions for the partonic off-shell 
matrix elements which were obtained in our previous paper [24]. 

Also we investigate here the prompt photon and associated muon production at Fermilab 
Tevatron. In the Standard Model (SM) these events come mainly due to Compton scattering 
process g + Q — > 7 + Q, where final state muon originates from the semileptonic decay of 
heavy (charm or bottom) quark Q [5]. It is important that studying of such processes 
possible can provide an information about new physics beyond the SM [28]. Therefore it is 
necessary to have a realistic estimation of associated 7 + /i production cross sections within 
the QCD. Such calculations in the /c T -factorization QCD approach are performed for the first 
time. In order to investigate the underlying dynamics in more detail, we study the azimuthal 
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correlations between the transverse momenta of produced prompt photon and muon. These 
quantities are sensitive to the different production mechanisms and also are powerful tests 
for the non-collinear evolution [29, 30]. 

The our paper is organized as follows. In Section 2 the KMR unintegrated parton densities 
in a proton are obtained and their properties are discussed. In particular, we compare the 
KMR gluon distributions with ones obtained from the full CCFM equation and within the 
framework of the Linked Dipole Chain model [31] (which is reformulation and generalization 
of the CCFM model). In Section 3 we present the basic formulas with a brief review of 
calculation steps. In Section 4 we present the numerical results of our calculations. Finally, 
in Section 5, we give some conclusions. 

2 The KMR unintegrated partons 

The Kimber-Martin-Ryskin approach [22] is the formalism to construct parton distri- 
butions f a (x, k^,yU 2 ) unintegrated over the parton transverse momenta from the known 
conventional parton distributions a(x,fj, 2 ), where a = xg or a = xq. This formalism is valid 
for a proton as well as photon and can embody both DGLAP and BFKL contributions. It 
also accounts for the angular ordering which comes from coherence effects in gluon emission. 
The key observation here is that the \x dependence of the unintegrated parton distributions 
f a (x, ky, (j, 2 ) enters at the last step of the evolution, and therefore single scale evolution 
equations (DGLAP or unified DGLAP-BFKL [32]) can be used up to this step. Also it was 
shown [22] that the unintegrated distributions obtained via unified DGLAP-BFKL evolution 
are rather similar to those based on the pure DGLAP equations. It is because the imposition 
of the angular ordering constraint is more important than including the BFKL effects. Based 
on this point, in our calculations we use much more simpler DGLAP equation up to the last 
evolution step. In this approximation, the unintegrated quark and gluon distributions are 
given [22] by 




where P a b(z) are the usual unregulated leading order DGLAP splitting functions, and q(x, fi 2 ) 
and g(x,fi 2 ) are the conventional quark and gluon densities. The theta functions which 
appear in (1) and (2) imply the angular-ordering constraint A = /i/(/i+ |k T |) specifically 
to the last evolution step to regulate the soft gluon singularities. For other evolution steps, 
the strong ordering in transverse momentum within the DGLAP equations automatically 
ensures angular ordering. It is important that parton distributions f a (x, k^, /x 2 ) extended 
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now into the k|. > /i 2 region. This fact is in the clear contrast with the usual DGLAP 
evolution 1 . 

The virtual (loop) contributions may be resummed to all orders by the quark and gluon 
Sudakov form factors 



\nT q (k 2 T ,n 2 ) = - J 



dp 2 T a s (pl) 



dzP qq (z), 



(3) 



A 1 

lnT 9 (k^ 2 ) = - j 



dprO! s (Pr) 



Pt 



2tt 



J- <-max 
71/ y dzPgg(z) + j dZZPgg(z) 



(4) 



where z max = 1 — <2 m in = /•*/(/•* + |Pt|)- The form factors T a (k^,/i 2 ) give the probability of 
evolving from a scale k 2 - to a scale /x 2 without parton emission. In according with (3) and 
(4) T a (kf,, /j 2 ) = 1 in the kf, > /j, 2 region. 

Note that such definition of the / a (x,k|,/i 2 ) is correct for k| > /i 2 , only where /i ~ 1 
GeV is the minimum scale for which DGLAP evolution of the collinear parton densities 
is valid. Everywhere in our numerical calculations we set the starting scale /i to be equal 
Ho = 1 GeV. Since the starting point of this derivation is the leading order DGLAP equations, 
the unintegrated parton distributions must satisfy the normalisation condition 



a(x, /j, 



fa(x, kj>, /i )d\s. T . 



(5) 



This relation will be exactly satisfied if we define [22] 



f a (x, k T , \i 



a(x,// 2 )T a (/i 2 ,/i 2 ). 



(6) 



Then, we have obtained the unintegrated parton distributions in a proton. In Figure 1 we 
show these densities f a (x, k 2 -, /i 2 ) at scale /i 2 = 100 GeV 2 as a function of x for different 
values of k 2 ,, namely k 2 , = 2 GeV 2 (a), k 2 T = 5 GeV 2 (b), k 2 T = 10 GeV 2 (c) and k 2 T = 
20 GeV 2 (d). The solid, dashed, short dashed, dotted, dash-dotted and short dash-dotted 
lines correspond to the unintegrated gluon (divided by factor 10), u+u, d+d, s, c and b quark 
distributions, respectively. We have used here the standard leading-order Gluck-Reya-Vogt 
(GRV) parametrizations [33] of the collinear quark and gluon densities a(x, (j 2 )- Note that 
unintegrated c and b quark distributions at k 2 ^ = 2 GeV 2 are very similar and cannot be 
really resolved in Figure 1. Also we have checked numerically that normalization condition 
(5) is correctly satisfied for all unintegrated parton distributions /„( 

It is interesting to compare the KMR-constructed unintegrated parton densities with the 
distributions obtained in other approaches. Recently the full CCFM equation in a proton 
was solved numerically using a Monte Carlo method, and new fits of the unintegrated gluon 
distributions (J2003 set 1 — 3) have been presented [34]. The input parameters were fitted 

1 We would like to note that cut-off A can be taken A = |1ct|/V also [25]. In this case the unintegrated 
parton distributions given by (1) — (2) vanish for > fi 2 in accordance with the DGLAP strong ordering 
in ky. 
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to describe the proton structure function F 2 (x, Q 2 ). These unintegrated gluon densities were 
used also in analysis of the forward jet production at HERA, charm and bottom production 
at Tevatron [34], and charm and J ftp production at LEP2 energies [35]. Also three different 
versions of the unintegrated gluon distributions obtained in the framework of the Linked 
Dipole Chain (LDC) model have been presented [36]. These gluon densities has been fitted 
to the inclusive F 2 data at HERA and already were used, in particular, in analysis [37] of 
the charm and beauty hadroproduction at Tevatron. In Figure 2 we plot the KMR (as a 
solid lines), the J2003 set 1 (as a dashed lines) unintegrated gluon distributions and so called 
standard version of the LDC unintegrated gluon distribution (as a dash-dotted lines) at scale 
ji 2 = 100 GeV 2 as a function of x for different values of gluon k^, namely k^ = 2 GeV 2 (a), 
\t 2 T = 10 GeV 2 (b), k 2 ; = 20 GeV 2 (c) and \t 2 T = 50 GeV 2 (d). One can see that behaviour of 
different unintegrated parton distributions in the small kf, region (which essentially drives 
the total cross sections) is very different. At the same time the differences between these 
distributions tends to be small when gluon transverse momentum k 2 - is large. Therefore the 
dependence of our predictions on the evolution scheme possible may be rather large, and 
further theoretical attempts are necessary to reduce this uncertainty. 



3 Calculation details 

3.1 The subprocesses under consideration 

The main contribution to the prompt photon production in proton-antiproton collisions at 
Tevatron and LHC colliders comes from quark-gluon and quark-antiquark induced partonic 
subprocesses: 

q {k l )+g{k 2 )^ 1 {f) + q{ P '), (7) 

q(kl) + q(k 2 )^'y(p^+g(p , ), (8) 

where the particles four-momenta are given in parentheses. Additionally, photons can be 
produced through the fragmentation of a partonic jet into a single photon carrying a large 
fraction z of the jet energy [7]. These processes are described in terms of quark-to-photon 
D q ^ 1 (z, ji 2 ) and gluon-to-photon D g ^(z, /i 2 ) fragmentation functions [8]. The main feature 
of the fragmentation contribution in leading order is fact that produced jet is balanced by a 
photon on the opposite side of the event and accompanied by collinear photon on the same 
side of the event. 

It is important that in order to reduce the huge background from the secondary photons 
produced by the decays of tt° and rj mesons the isolation criterion is introduced in the 
experimental analyses. This criterion is the following. A photon is isolated if the amount of 
hadronic transverse energy E^ ad , deposited inside a cone with aperture R centered around 
the photon direction in the pseudo-rapidity and azimuthal angle plane, is smaller than some 
value £f? ax : 

Tjihad ^ Tjimax 

T ' (9) 

The both D0 and CDF collaborations take R ~ 0.4 and Z?™ ax ~ 1 GeV in the experiment [1- 
5] . Isolation not only reduces the background but also significantly reduces the fragmentation 
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components. It was shown [38] that after applying the isolation cut (9) the contribution 
from the fragmentation subprocesses is about 10% of the total cross section. Since the 
dependence of our results on the non-collinear parton evolution scheme may be rather large 
(as it was demonstrated in Section 2), in our further analysis we will neglect the relative 
small fragmentation contribution and consider only the q + g 7 + q and q + q — > 7 + g 
subprocesses. We would like to note that photon produced in these processes is isolated from 
the quark or gluon jet by requiring a non-zero transverse momentum of a photon or jet in 
the pp center-of-mass frame. 



3.2 Kinematics 

Let pi and P2 be the four-momenta of the incoming protons. The initial off-shell par- 
tons have the four-momenta k\ and k 2 . In our analysis below we will use the Sudakov 
decomposition, which has the following form: 

p 1 = aipi + (3ip 2 + Pt, p' = OL 2 p\ + P2P2 + p' T , (7) 
h = x x p\ + fcir, k 2 = X2P2 + k 2T , 

where k±T, k 2 r, Pt an d Pt are the transverse four- momenta of the corresponding particles. 
It is important that k\ T = — k\ T 7^ and k\ T = —k\ T 7^ 0. In the pp center-of-mass frame 
we can write 

01 = ^/2(1,0,0,1), ^2 = ^/2(1,0,0,-1), (8) 

where s — (pi + p 2 ) 2 is the total energy of the process under consideration and we neglect 
the masses of the incoming protons. The Sudakov variables are expressed as follows: 

«i = ^7=exp(y 7 ), a 2 = ^pexp(y'), 
V s V s 



^1 m > (9) 

/3i = exp(-y 7 ), p 2 = —L exp(-y'), 



e t t i\ a m 
V s V s 
where E T and m' T are the transverse enegy and transverse mass of produced photon and 
outgoing parton, respectively, and y 1 and y' are their rapidities (in the pp center-of-mass 
frame). The photon pseudo-rapidity rp is defined as rf 1 = — lntan(# 7 /2), where 9 1 is the 
polar angle of the prompt photon with respect to the proton beam. From the conservation 
laws we can easily obtain the following conditions: 

x 1 =a 1 + a 2 , x 2 = (3 1 +(3 2 , k iT + k 2T = Pt + Pt- ( 10 ) 

The scaling variable xt = 2E T /^ is also often used in analysis of the prompt photon 
production. 



3.3 Cross section for prompt photon production 

The total cross section for prompt photon hadroproduction at high energies in the kr- 
factorization QCD approach can be written as 



77 J Xi Z7T 

a,b L 



X 



/ — fb{x 2 , \4 T , fi 2 )dkl T — -dcr(ab -> 7c), 
J x 2 2ir 



(11) 
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where a. . . c = q and/or g, a(ab — > 7c) is the cross section of the photon production in 
the corresponding parton-parton interaction. Here initial partons a and b have longitudinal 
momentum fractions X\ and x 2 , non-zero transverse momenta k 1T and k 2T and azimuthal 
angles 0i and (f> 2 . From the expression (11) we can easily obtain the formula: 

, (p + ^ 7 + X) = E/^ ? l^lV-7c)x 

a ,b v ; ^ 



xf a (xi,kl T ,fi 2 )f b (x 2 ,kl T ,ii 2 )dkl T dkl T dE^dy y dy c 



2tt 2tt 2tt ' 



where |7W | 2 (a6 — > 7c) is the hard partonic squared off-mass shell matrix element which 
depends on the transverse momenta k\ T and k 2 T , y c is the rapidity of the parton c in the pp 
center-of-mass frame, 0i, (fi 2 and 7 are the azimuthal angles of the initial partons a, b and 
produced prompt photon, respectively. The analytic expression for the \M.\ 2 (ab — > 7c) was 
obtained in our previous paper [24] . We would like to note that if we average the expression 
(12) over k\T and k 2 x and take the limit k\ T — > and k 2T — > 0, then we obtain well-known 
expression for the prompt photon hadroproduction in leading-order (LO) perturbative QCD. 

The multidimensional integration in (12) has been performed by means of the Monte 
Carlo technique, using the routine VEGAS [39]. The full C++ code is available from the 
authors on request 2 . 



4 Numerical results 

We now are in a position to present our numerical results. First we describe our theoret- 
ical input and the kinematical conditions. After we fixed the unintegrated parton distribu- 
tions in a proton f a (x, k\, /z 2 ), the cross section (12) depends on the energy scale \i. As it 
often done [38] for prompt photon production, we choose the factorization and renormaliza- 
tion scales to be equal fip = fiR = fi = ^Ej,. In order to estimate the theoretical uncertainties 
of our calculations we will vary the scale parameter £ between 1/2 and 2 about the default 
value £ = 1. Also we use LO formula for the strong coupling constant a s (fi 2 ) with n/ = 3 
active (massless) quark flavours and Aq C d = 232 MeV, such that a s (M|) = 0.1169. In our 
analysis we not neglect the charm and bottom quark masses and set them to be m c = 1.5 
GeV and m h = 4.75 GeV, respectively. 

4.1 Inclusive prompt photon production at Tevatron 

Experimental data [1-4] for the inclusive prompt photon hadroproduction p + p — > 7 + X 
come from both the D0 and CDF collaborations. The D0 [1, 2] data were obtained in the 
central and forward pseudo-rapidity regions for two different center-of-mass energies, namely 
y/s = 630 GeV and y/s = 1800 GeV. The central pseudo-rapidity region is defined by the 
requirement \r]"'\ < 0.9, and the forward one is defined by 1.6 < \r]~'\ < 2.5. The more 
recent CDF data [3] refer to the same central kinematical region |r] 7 | < 0.9 for both beam 
energies y/s = 630 GeV and y/s = 1800 GeV. Also very recently the CDF collaboration has 

2 lipat ov@ theory, sinp .msu.ru 
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presented a new measurement [4] of the prompt photon cross section at ^/s = 1800 GeV. 
This measurement is based on events where the photon converts into an electron-positron 
pair in the material of inner detector, resulting in a two track event signature ("conversion" 
data). These data refer only to the central kinematical region. Actually, there are two 
different datasets, which were used in the CDF measurement with conversions, namely 8 GeV 
electron data and 23 GeV photon data 3 . In all these measurements the double differential 
cross sections da / dE^drp as a function of the transverse energy Ej, are determined. 

The results of our calculations are shown in Figs. 3 — 9. So, Figs. 3 and 4 confront the 
cross sections da/dE^drf calculated at y/s = 630 GeV in different kinematical regions with 
the D0 [2] and CDF [3] data. The solid lines are obtained by fixing both the factorization 
and renormalization scales at the default value fx = Ej, whereas upper and lower dashed 
lines correspond to the fx = E^/2 and fx = 2Ej, scales, respectively. One can see that our 
predictions agree with the experimental data within the scale uncertainties. However, the 
results of calculation with the default scale tend to underestimate the data in the central 
kinematical region and agree with the D0 data in the forward rp region. The collinear NLO 
QCD calculations [38] give the similar description of the data: the results of measurement are 
higher than the NLO prediction at low in the central rp range but agree at all E? in the 
forward pseudo-rapidity region. Then, one can see that the scale dependence of our results is 
rather large: the variation in fx as it was described above changes the cross sections by about 
20 — 30%. The theoretical uncertainties of the collinear NLO calculations are similar (about 
20%) [2]. However, one should keep in mind that additional dependence of our results on 
the evolution scheme may be also rather large (as it was discussed in Section 2), and overall 
agreement with the experimental data can be improved when unintegrated quark and gluon 
distributions in a proton will be studied more detail. At the same time the use of different 
sets of the parton distributions in NLO calculations changes the cross sections by less than 
5% [2, 3]. 

The double differential cross sections da/dE^drf compared with the experimental data at 
\fs = 1800 GeV in different pseudo-rapidity regions are shown in Figs. 5 — 7. All curves here 
are the same as in Fig. 3. We find that our predictions agree well with the D0 [1] and CDF [3, 
4] data both in normalization and shape. There are only rather small overestimation of the 
data at low Ej. values in Figs. 6 and 7. Again, the scale dependence of our calculations 
is about 20 — 30%. The theoretical uncertainties of the collinear NLO predictions are 
much smaller, about 6% [1]. Note that now NLO calculations agree with the data more 
qualitatively. So, the shape of the measured cross sections is generally steeper than that of 
the NLO predictions. It was shown [3, 4] that this shape difference is difficult to explain 
simply by changing the renormalization/factorization scales in the collinear calculation, or 
the set of parton distribution functions. 

Also the disagreement between data and NLO calculations is visible [2, 3] in the ratio of 
the cross sections at different energies. This quantity is known as a very informative subject 
of investigations and provides a precise test of the QCD calculations. It is because many 
factors which affect the absolute normalization, as well as many theoretical and experimental 
uncertainties partially or completely cancel out [2, 3]. In particular, the cross section ratio 
provides a direct probe of the matrix elements of the hard partonic subprocesses since the 

3 See Ref. [4] for more details. 
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theoretical uncertainties due to the quark and gluon distributions are reduced. 

So, the D0 collaboration has published the results of measurement [2] for the ratio of 
630 GeV and 1800 GeV dimensionless cross sections a D as a function of scaling variable x T - 
The measured cross section cr D averaged over azimuth is defined as 



The ratio <td(630 GeV)/(T£)(1800 GeV) compared with the D0 experimental data [2] in dif- 
ferent pseudo-rapidity rp regions is shown in Figs. 8 and 9. The solid lines represent the 
/c T -factorization predictions at default scale /i. For comparison we show also the results of 
the collinear leading-order (LO) QCD calculations with the GRV parton densities [33] of 
a proton (as a dashed lines). Note that when we perform the LO QCD calculations we 
take into account the partonic subprocesses (7) and (8) and neglect the small fragmenta- 
tion contributions, as it was done in the /cy-factorization case. It is clear that although the 
experimental points have large errors they tend to support the /c T -factorization predictions. 
We would like to point out again that now sensitivity of our results to the non-collinear 
evolution scheme is minimized. In the collinear approach, the NLO corrections improve the 
description of the data and then sum of LO and NLO contributions practically coincides with 
our results at x T > 0.05 [2]. This fact is clear indicates that the main part of the collinear 
high-order corrections is already included at leading-order level in the /c T -factorization for- 
malism 4 . Nevertheless, the experimental data at the lowest xt are systematically higher [2] 
than NLO QCD predictions in both central and forward pseudo-rapidity regions, and it was 
claimed [3] that ratio of the two cross sections as a function of the xt is difficult to reconcile 
with the NLO QCD calculations. 

4.2 Associated prompt photon and muon production at Tevatron 

Now we investigate the prompt photon and associated muon production at Tevatron. 
These events are assumed to come from associated prompt photon and heavy (charm or 
bottom) quark production with the heavy quark decaying into a muon [5]. To calculate 
muon production from heavy quarks we first convert charm (bottom) quarks into D (B) 
hadrons using the Peterson fragmentation function [40] and then simulate their semileptonic 
decay according to Standard Model. Our default set of the fragmentation parameters and 
branching rations is e c = 0.06, f(c -> //) = 10.2% and e b = 0.06, f(b -> fi) = 10.8%. 
Of course, muon transverse momenta spectra are sensitive to the fragmentation functions. 
However, this dependence is expected to be small as compared with the scale uncertainties 
and the uncertainties connected with unintegrated parton densities. On the other hand, the 
variations of e c and e& do not affect on the azimuthal angle distribution (which is one of the 
main subject of our study) because a fragmentation does not change the direction of the 
quark or hadron momentum. 

The experimental data for the 7 + /i cross section at Tevatron come from CDF collabora- 
tion [5]. The differential cross section da/dp? at yHi = 1800 GeV was obtained. The photon 
pseudo-rapidity is required to be within \rp\ < 0.9 region whereas the muon transverse 
momentum and pseudo-rapidity are required to be > 4 GeV and |r/ M | < 1.0. 

4 See Refs. [15, 16] for more details. 
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The transverse momentum distribution da/dpj, in comparison to experimental data [5] 
is shown in Fig. 10. All curves here are the same as in Fig. 3. One can see that the 
shape of this distribution is well described by our calculations. However, the theoretical 
results slightly overestimate the data in absolute normalization. This fact can be connected 
with the GRV [33] parametrization of charm and bottom collinear densities of a proton 
(which are one of the basic ingredients of our derivation). Note, however, that in general 
the experimental points still lie within scale uncertainties (about 30%) of our calculations. 
It is important also that our predictions practically coincide with the results of collinear 
NLO QCD calculations [41], which are much larger than LO ones [5]. This fact clearly 
demonstrates again that the main part of the standard high-order corrections is already 
effectively included in the /^-factorization approach. 

Further understanding of the process dynamics and in particular of the high-order effects 
may be obtained from the angular correlation between the transverse momenta of the final 
state particles. It was shown [30] that investigation of these correlations is a powerful test for 
the non-collinear parton evolution dynamics. This is because such quantities are sensitive 
to relative contributions of different production mechanisms to the total cross section. So, 
in the collinear LO approximation, the prompt photon and heavy quark Q are produced 
back-to-back. Therefore distribution over the azimuthal angle difference A0 7< ^ must be 
simply a delta function at A0 7<3 = n. Taking into account the non-vanishing initial parton 
transverse momenta k 1T and k 2 r leads to the violation of this back-to-back kinematics in 
the /cT-factorization approach. 

The differential cross section da/dA^ 1 ^ calculated at p^ > 4 GeV, \r]^\ < 1.0 and < 
0.9 is shown in Fig. 11. The solid and both dashed lines here are the same as in Fig. 3. The 
LO QCD contribution also shown (as a dash-dotted line). One can see a clear difference 
in shape between /c^-factorization results and collinear LO QCD ones. As it was expected, 
the LO QCD predicts only a peak at A0 W = n. The small broadering of A0 7M distribution 
illustrates the effect of Q — > \x decays. Unfortunately, the predictions of the NLO QCD 
for this distribution are still unknown. The direct comparison between NLO calculations 
and our results should give a number of interesting insighs. In particular, it can provide 
an important information about high-order terms which are missed in the leading-order kr- 
factorization approach. In any case, the future theoretical and experimental study of such 
processes will be important check of non-collinear parton evolution. 

4.3 Inclusive prompt photon production at LHC 

We can conclude from Figs. 3 - - 11 that our calculations in general agree well with 
experimental data [1-5] taken by the D0 and CDF collaborations at Tevatron. Based on 
this point, now we can try to extrapolate our theoretical predictions to LHC energies. We 
perform the calculation for both central and forward pseudo-rapidities rp . As a representative 
example, we will define the central and forward kinematical regions by the requirements 
\r] 1 \ < 2.5 and 2.5 < \r] 1 \ < 4, respectively. 

The transverse energy Ej, distributions of the inclusive prompt photon production in 
different if ranges at ^[s = 14 TeV are shown in Figs. 12 and 13. All curves here are the 
same as in Fig. 3. One can see that variation in scale \i changes the estimated cross sections by 
about 20 — 30%. However, as it was already discussed above, there are additional theoretical 
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uncertainties due to the non-collinear parton evolution, and these uncertainties are not well 
studied up to this time. Also the extrapolation of the available parton distribution to the 
region of lower x is a special problem at the LHC energies. In particular, one of the problem 
is connected with the correct treatment of saturation effects in small x region 5 . Therefore 
much more work needs to be done until these uncertainties will be reduced. 

5 Conclusions 

We present calculations of the prompt photon hadroproduction at high energies in the 
/cT-factorization QCD approach. In order to obtain the unintegrated quark and gluon dis- 
tributions in a proton we have used the Kimber-Martin-Ryskin prescription. We have in- 
vestigated both inclusive prompt photon and associated with muon production rates. The 
associated 7+// events come primarily due to the Compton scattering process g + Q — * 7 + Q, 
with the final state heavy (charm or bottom) quark Q producing a muon. The calculations 
of such cross sections in the /^-factorization approach were performed for the first time. 

We have found that our predictions for the inclusive prompt photon production agree well 
with experimental data taken by the D0 and CDF collaborations at Tevatron in both central 
and forward pseudo-rapidity regions. It is very important that perfect agreement was found 
also in the ratio of two cross sections calculated at y/s = 630 GeV and y/s = 1800 GeV. This 
ratio provides a direct probe of the QCD matrix elements since in this case the theoretical 
uncertainties connected to the parton distributions are significantly reduced. We have also 
demonstrated that the main part of the standard high-order corrections is already included 
in the /c^-factorization formalism at LO level. Additionally, we present our predictions for 
the inclusive prompt photon production at LHC. 

At the same time our results for associated 7 + /i production tend to overestimate the 
CDF data but still agree with data within the scale uncertainties. We have demonstrated 
also the significant deviation from back-to-back kinematics in prompt photon and associated 
muon production. We can expect that further theoretical and experimental study of such 
processes will give an important information about non-collinear parton evolution dynamics. 

In order to investigate the theoretical uncertainties of our results we have studied the 
sensitivity of our predictions to the renormalization and factorization scales. We have found 
that this dependence is about 20 — 30% in wide center-of-mass energy range. There are, 
of course, also additional uncertainties due to unintegrated parton distributions in a proton. 
Therefore much more work needs to be done before these uncertainties will be reduced. 
Finally, in our analysis we neglect the contribution from the fragmentation processes. We 
plan to investigate these problems in more detail in the forthcoming publications. 
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Figure 1: The unintegrated parton distributions f a ( ) at scale \i 2 = 100 GeV 2 as 

a function of x for different values of k|>, namely k|> = 2 GeV 2 (a), k 2 ^ = 5 GeV 2 (b), 
k 2 ; = 10 GeV 2 (c) and k| = 20 GeV 2 (d). The solid, dashed, short dashed, dotted, dash- 
dotted and short dash-dotted lines correspond to the unintegrated gluon (divided by factor 
10), u + u, d + d, s, c and b quark distributions, respectively. 
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Figure 2: The unitegrated gluon distributions f g ( ) at scale /z 2 = 100 GeV 2 as a 

function of x for different values of k 2 ^, namely k 2 ^ = 2 GeV 2 (a), k^ = 10 GeV 2 (b), k^ = 
20 GeV 2 (c) and k^ = 50 GeV 2 (d). The solid, dashed and dash-dotted lines correspond 
to the KMR, J2003 set 1 and standard version of the LDC unintegrated gluon densities, 
respectively. 
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Figure 3: The double differential cross section da/dE^drf for inclusive prompt photon 
hadroproduction at \rp\ < 0.9 and y/s = 630 GeV. The solid line corresponds to the de- 
fault scale jj, = Ej,, whereas upper and lower dashed lines correspond to the \i = and 
ji = IE scales, respectively. The experimental data are from D0 [2] and CDF [3]. 
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Figure 4: The double differential cross section dajdE^drp for inclusive prompt photon 
hadroproduction at 1.6 < |r? 7 | < 2.5 and y/s = 630 GeV. All curves are the same as in 
Figure 3. The experimental data are from D0 [2]. 
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Figure 5: The double differential cross section da/dE^drf for inclusive prompt photon 
hadropro duct ion at \rf'\ < 0.9 and y/s = 1800 GeV. All curves are the same as in Fig- 
ure 3. The experimental data are from D0 [1] and CDF [3]. 
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Figure 6: The double differential cross section dajdE^drp for inclusive prompt photon 
hadroproduction at 1.6 < \rf 1 \ < 2.5 and y/s = 1800 GeV. All curves are the same as 
in Figure 3. The experimental data are from D0 [1]. 
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Figure 7: The double differential cross section da/dB^drf for inclusive prompt photon 
hadroproduction at |r? 7 | < 0.9 and y/s = 1800 GeV. All curves are the same as in Fig- 
ure 3. The experimental data are from CDF [4]. 
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Figure 8: The ratio of the dimensionless cross sections a D (630 Ge V)/ctd (1800 GeV) as a func- 
tion of scaling variable xt at \r] 1 \ < 0.9. The solid line was obtained in the /c T -factorization 
approach whereas dashed line corresponds to the collinear leading-order QCD calculations 
with the GRV parton densities of the proton. The renormalization and factorization scales 
are \x = Ej.. The experimental data are from D0 [2]. 
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Figure 9: The ratio of the diinensionless cross sections o D (630GeV)/o D (1800GeV) as 
function of scaling variable X? at 1.6 < \if'\ < 2.5. All curves are the same as in Figure 
The experimental data are from D0 [2] . 
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Figure 10: The differential cross section da/dpj, for associated prompt photon and muon 
hadroproduction at |r? 7 | < 0.9, |^| < 1.0, > 4 GeV and ^/s = 1800 GeV. All curves are 
the same as in Figure 3. The experimental data are from CDF [5]. 
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Figure 11: Azimuthal correlations in associated prompt photon and muon hadroproduction 
at \rp\ < 0.9, |^| < 1.0, Vt > 4 GeV and = 1800 GeV - The solid line corresponds to 
the default scale fi = Ej>, whereas upper and lower dashed lines correspond to the fi = 
and \i = IEj. scales, respectively. The dash-dotted line correspond to the sole LO QCD 
calculations at /i = E^. 
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Figure 12: The differential cross section da /dEj, for inclusive prompt photon hadropro emo- 
tion at \r] 1 \ < 2.5 and y/s = 14 TeV. All curves are the same as in Figure 3. 
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Figure 13: The differential cross section da /dEj, for inclusive prompt photon hadropro emo- 
tion at 2.5 < \r] 1 \ < 4 and y/s = 14 TeV. All curves are the same as in Figure 3. 
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